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It has recently been shown that specific non-perturbative effects may lead to an explosive decay 
of flat direction condensates in supersymmetric theories. We confirm explicitly the efficiency of this 
process with lattice simulations: after few rotations of the condensates in their complex plane, most 
of their energy is quickly converted into inhomogeneous fluctuations. We then point out that this 
generates a gravitational wave background which depends on the inflaton sector and falls in the 
Hz-kHz frequency range today. We compute the resulting spectrum and study how it depends on 
the parameters. We show that these gravity waves can be observable by upcoming experiments like 
Advanced LIGO and depend crucially on (i) the initial VEV of flat directions when they start to 
oscillate, (ii) their soft SUSY-breaking mass and (iii) the reheat temperature of the universe. This 
signal could open a new observational window on inflation and low-energy supersymmetry. 
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Gravity wave (GW) experiments could provide unique 
informations about high-energy phenomena in the early 
universe, such as inflation, preheating after inflation, cos- 
mic strings and first order phase transitions In par- 
ticular, cosmological GW backgrounds will be a target 
for several high-sensitivity interferometric experiments, 
planned or proposed, in the frequency range between 
10 -5 and 10 3 Hz. In this letter, we want to consider a 
new cosmological source of GW, which emerges naturally 
in the framework of high-energy physics. 

Supersymmetric theories typically involve many flat 
directions (see @ for the MSSM), i.e. directions in 
field space where the renormalizable part of the scalar 
potential is exactly flat in the limit of unbroken su- 
persymmetry. The flatness is lifted by soft SUSY- 
breaking terms, non-renormalizable terms and SUSY- 
breaking terms from the finite energy density in the early 
universe [3j. Scalar field condensates may develop large 
VEVs along these directions during inflation and have 
several interesting consequences in cosmology Q ■ A ma- 
jor example is the Affleck-Dine mechanism for baryoge- 
nesis @ . Condensates with large VEVs start to oscillate 
when the Hubble rate H becomes of the order of their 
soft mass m ~ TeV (lij . with an initial amplitude that 
is model-dependent and could be as high as the Planck 
scale [3, 8]. It has recently been shown [j| that 
non-perturbative resonant effects may lead to an explo- 
sive decay of these coherent oscillations. We will show 
that this generates GW that can be observable by up- 
coming experiments fl2j . 

The resonant decay of flat directions shares many simil- 
itudes with preheating after inflation, when the inflaton 
decays into large, non-thermal fluctuations of itself and 
other bosonic fields [l4[ ■ GW from preheating have been 
intensively studied recently [HI, Ha. 1171 Il8j . Ref. [15[ de- 
veloped theoretical and numerical methods to calculate 
GW production from dynamical scalar fields in an ex- 
panding universe. We will apply and extend these meth- 
ods to the decay of flat directions. 

The resulting GW, however, will have different prop- 
erties. GW from preheating have been studied so far in 



two main classes of models: parametric resonance after 
chaotic inflation and tachyonic preheating after hybrid 
inflation. In the first case, the typical GW frequency 
today is of the order of / ~ 10 7 Hz, which is too high 
to be observable. Hybrid inflation models may occur at 
lower energy scales, so it was initially hoped that GW 
from tachyonic preheating may fall into an observable 
range. However, as was shown in [T^|. this generically 
requires very small coupling constants [191 ]. By contrast, 
we will see that GW from the non-perturbative decay 
of flat directions fall naturally in the Hz-kHz frequency 
range today, where high-sensitivity interferometric ex- 
periments can operate. Another difference - which will 
be crucial for the GW amplitude - is that, whereas the 
inflaton dominates the energy density before it decays, 
flat directions are usually subdominant (at least for sub- 
Planckian VEVs and when their non-perturbative decay 
is sufficiently fast). During their decay, the universe can 
then be matter-dominated if the inflaton is still oscil- 
lating around the (quadratic) minimum of its potential, 
or radiation-dominated if it decays earlier and its decay 
products thermalize when H > m. Each case leads to 
different GW spectra, so GW from the non-perturbative 
decay of flat directions carry also informations on the 
inflaton sector. 

To illustrate why it may be difficult for GW from 
preheating to be observable, consider a simple model, 
V = m 2 (j) 2 + g 2 <f) 2 x 2 , W1 th two real scalar fields 4> and x- 
In this model, %-particles are produced due to the non- 
adiabatic evolution of their frequency as the condensate 
6 oscillates around the minimum UJ. For to - TeV, 
this could a priori describe the decay of a flat direction 
4>. We will denote by an index "i" the time U when the 
flat direction starts to oscillate (Hi = to), by <&i its ini- 
tial amplitude at that time and we normalize the scale 
factor to at = 1. The spectra of the quanta amplified 
by preheating are usually strongly peaked around some 
typical (comoving) momentum fc*, which is inherited by 
the GW spectru m flM [l6l Ha |. For the model above, 
fc* ~ y/gmQid 1 / 4 [Tj|. Assuming for simplicity that the 
universe is already radiation-dominated during the de- 
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cay, this leads to the estimate /* ~ \fg^]Mp\ 10 11 Hz 
for the peak frequency of the resulting GW today. Tak- 
ing for instance $ = 1CP 2 Mpi, /* < 10 3 Hz requires an 
extremely small coupling constant g 2 < 10 -28 . If <&i is 
smaller, the GW amplitude is too low to be observable. 
Indeed, since the energy density in the flat direction p^t 
is subdominant, only a small fraction pnat/ptot of the to- 
tal energy density is available for the production of GW. 
We will see that the fraction of energy density in GW 
varies as (pa a t / Ptot) 2 , which can be very small. 

However, flat directions are complex fields and non- 
renormalizable terms usually generate a velocity for their 
phase when H ~ to 0, Q. The subsequent dynamics 
is dominated by the m 2 \4>\ 2 term, resulting in a con- 
stant, non-zero phase velocity and an elliptical motion for 
a 3 / 2 in its complex plane. This leads to time-dependent 
(and quasi-periodic) mixings between different excita- 
tions around flat directions, which may result in a very 
efficient mechanism of particle production 0, [13] , signifi- 
cantly different from the model considered above. In par- 
ticular, these particles have a typical momentum fc* ~ to, 
of the order of the Hubble scale. Formally, this amounts 
to replace y/g<f>i by ^fm in the previous model. This gives 
a peak frequency of GW that is independent of <I>i and 
of order 10 3 Hz today for m ~ TeV. 

This encouraging estimate should be complemented by 
detailed computations. The resonant decay of flat direc- 
tions has been studied so far only in the linear regime, ne- 
glecting any backreaction of the quanta produced. Back- 
reaction is eventually responsible for most of the decay 
and we have to use lattice simulations to fully take it into 
account. We will simulate the model d, |2(j 

v = to 2 h 2 + to 2 \x? + 9 ^ {4>x* + <t>*x) (i) 

where * denotes the complex conjugate, <f> the flat direc- 
tion field and x its decay products. This potential in- 
cludes the relevant terms for excitations around MSSM 
flat directions (where the interaction comes from the D- 
term) but at least two flat directions have to be consid- 
ered when the symmetry broke n by their VEV is gauged 
[l|. This case was studied in [nj, [ll| and found to be 
largely similar to the model (Q]), so here we focus on (fT|) 
for simplicity. We perform lattice simulations with the 
parallel code ClusterEasy [2l| and we compute the re- 
sulting GW spectrum with the method of [lq |. extended 
here to take into account the expansion of the universe. 

Consider first the non-perturbative decay itself. The 
background evolution of the condensate is given by ip = 
(cosTOi + ie sinTOi) <I> l /\/2a 3 . The " ellipticity" e < 1 is 
model-dependent but typically of order e ~ 0.1 Q. An 
analytical study of particle production in this background 
will appear elsewhere. Here we just present the results. 
It is convenient to introduce the parameter q = qi/a 3 
with qi — g 2 Q 2 /m 2 . Note that qi can be as high as 
10 30 for <i>i close to the Planck scale. In the realistic 
case qe 4 » 1 [22), particle production occurs because 
the physical eigenstates vary with time Fig. [T] shows 



the spectra of the %-modes in a Minkowski background 
when backreaction is still negligible. The resonance is 
characterized by large instability bands separated by nar- 
row stability bands, contrary to parametric resonance but 
similarly to the tachyonic resonance studied in the sec- 
ond ref. of [HI]. When the trajectory of the condensate 
in its complex plane is circular, e = 1, there is a tachyonic 
amplification of the modes with k 2 + to 2 < m 2 9] . When 
e < 1 , the amplification is more efficient and extra insta- 
bility bands appear in the UV, roughly up to k < m/e. 
The IR band has the highest amplitude, except if m x is 
very close to m. The peak of the spectrum is located at 
k% = kz/m ~ 0.1. In an expanding universe, the phys- 
ical momenta k/a "move" in the resonance pattern and 
the peak is shifted to a comoving momentum *~ 1. 
When the %-modes have been sufficiently amplified, they 
backreact on the condensate and convert its energy into 
inhomogeneous fluctuations of both fields. This will be 
responsible for most of the GW production. We show the 
evolution of the mean of c?l 2 <\> in its complex plane in 
Fig. [2j The field first follows its elliptical trajectory, but 
the coherent motion is then quickly destroyed by backre- 
action. This typically happens after only 1 to 5 rotations 
of the condensate (the higher g 2 , the sooner it occurs). 
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FIG. 1: Spectra of the x modes without expansion of the 
universe when backreaction is still negligible. They are in- 
dependent of q as far as qe 4 ^> 1. For qe 4 1, particle 
production is much less efficient. 
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FIG. 2: Evolution of the mean of the real and imaginary parts 
of a 3 / 2 <j) in a matter-dominated background with e = 0.4, 
g 2 = 0.01, m = 1 TeV and $i = 10~ 12 Mpi. 

We now turn to the actual computation of GW. Their 
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evolution from the time of production (that we will de- 
note with an index "p") up to now depends on the back- 
ground equation of state, see e.g. [ill ] for details. In our 
case, their present-day frequency and energy density per 
logarithmic frequency interval [l| read 




where k is the comoving wave-number and pi is the 
total energy density at t t . If the universe is matter- 
dominated (MD) during the flat direction decay and be- 
comes radiation-dominated (RD) later, when a = cird, 
then GW are further redshifted and diluted by the fac- 
tors in a,i/aRD above. These factors are absent if the 
universe is already RD at t j [24[ ■ 

We can estimate how the GW spectrum depends on 
the parameters from the characteristic physical length 
i?* = a/k* of the scalar field inhomogeneities fl5l. HH [lq| . 
but we have to take into account that the flat direction 
energy density pa a t is usually subdominant. We then 
estimate the fraction of energy density in GW at the time 
of production as [H[ p gw /ptot ~ 0.1 (R*H) 2 (p flat / p to t) 2 ■ 
Using pfl a ti = to 2 $ 2 , fc* = k*to and Hi — m, Eqs.© give 

^--v^fe) 1/45x102112 ,3) 

~ V (ws (£) 

for the frequency and amplitude of the GW spectrum's 
peak. Note that the amplitude is very sensitive to <I>i. 
Note also that Eq.Q cannot be extrapolated to arbi- 
trary large $ j . If a flat direction acquires an initial VEV 
of order Mpi/3 or larger, it dominates the total energy 
density before decaying and its oscillations are delayed 
until $ l ~ Mpi/3. Thus Eq.© with $ t = M Pl /3 can be 
considered as an upper limit on the GW amplitude, i.e. 
h 2 Q* < 1CT 6 . The same bound follows from [H with 
R^^l/H and p flat <; p to t UM- 

We confirmed the predictions ([3l [4]) with intensive lat- 
tice simulations. Two different scales have to be kept 
under control numerically: the mass and the VEV of the 
condensate, whose ratio is measured by In our sim- 
ulations, to defines the natural unit of time while 1/x/Qi 
defines the time step required to accurately follow the 
dynamics. Clearly, we cannot simulate values of qi as 
high as 10 30 . However, we saw above that the linear 
stage of particle production is independent of g,-, see also 
Q. Indeed, the GW spectra that we obtained for dif- 
ferent qi (by varying g 2 ) were nearly identical, as far as 
qe 4 ^> 1 was satisfied. We were able to simulate val- 
ues of qi ranging from 10 4 to 10 8 , which allowed us to 
take e ranging from 0.2 to 1. Smaller values of e lead 
to faster particle production, which could increase the 
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FIG. 3: Lattice results for the GW spectra today, normalized 
as /i 2 f2 gw x Mp 1 /(10$i) 4 . Here we did not include the m/aRD 
factors in ©. Top-left: RD background, $; = 10~ 13 M P i, 
m = 122 GeV, m x = 0, qi = 10 6 , e = 0.5 (blue). Top-right: 
RD background, $i = 10 _1 M P i, m = 1.22 TeV, m x = 0, (i) 
qi = 10 6 , e = 0.5 (black) and (ii) q, = 10 8 , e = 0.2 (green). 
Bottom-left: MD background, = 10" 13 M P1 , m = 122 GeV, 
qi = 10 6 , (i) m x — 0, e — 0.5 (blue) and (ii) m x — m/5, 
e = 1 (black). Bottom-right: MD background, $i = 10 _1 M P i, 
m = 1.22 TeV, m x = 0,q t = 10 6 , e = 0.5 (red). 

GW amplitude. Increasing the mass of x h as th e °P" 
posite effect and makes the UV modes relatively more 
important (which makes this case more difficult to simu- 
late). However, this becomes negligible when m x differs 
significantly from to. This is probably more natural in 
models with multiple flat directions [lOj . We varied g 2 
from 1CT 24 to 0.5, from 10- 13 M P1 to 10 _1 M P1 and 
to from 100 GeV to 10 TeV, finding excellent agreement 
with Eqs. (OHJl, see Fig. [3] Note that ~ 1 is slightly 
smaller in a RD background than in a MD one, leading 
to GW with higher amplitude and smaller frequency at 
the time of production. 

In Fig.21 we compare our results with the sensitivity of 
interferometric experiments. We consider the soft masses 
to = 100 GeV and m = 10 TeV, different initial VEVs 
<I>i and different effective temperatures Tr oc p 1 ^ when 
the universe becomes radiation-dominated. If the infla- 
ton decays slowly (perturbatively) and its decay products 
quickly thcrmalizc, then Tr may correspond to the usual 
reheat temperature of the universe. In general, however, 
the universe may become radiation-dominated long be- 
fore full thermal equilibrium is established, see e.g. [23| . 
In particular, this means that the standard bound from 
the thermal over-production of gravitinos does not nec- 
essarily apply on Tr. For Tr J> 0.2y/mMp\, the universe 
is RD during the flat direction decay and the GW am- 
plitude is independent of m. If Tr is lower, GW are fur- 
ther diluted and redshifted by the factors in ai/aRu in 
Eqs. ([2]), which depend on to. Assuming a fast transition 
from MD to RD, we have ai/a RD - (5Tr/ 'VmMpjjf/ 3 in 
that case. We show the results for Tr — 10 9 GeV [27| and 
10 7 GeV. We see that these GW can be observable for a 
reasonable range of soft masses and Tr, but this gener- 
ically requires that some flat directions acquire a large 
enough initial VEV. A detection would be easier for low 
to or high Tr, but we should have at least <&i > 10 17 
GeV. 
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FIG. 4: Sensitivity sketch of Advanced LIGO and Einstein 
Telescope []| compared to GW spectra produced in the model 
|T]). In each case, the spectrum at the left is for m — 100 GeV 
and the one at the right for m = 10 TeV. 

The initial amplitude of flat direction oscillations de- 
pends on the lowest order of the non-renormalizable 
terms that lift the flatness of each direction [3|. A 
superpotential term like 4> n /M n ~ 3 (with M the cutoff 
scale) leads to $ 4 - (mF- 3 ) 1 /!"^) [3] w hen H ~ m. 
In the MSSM, n = 9 is the minimal value allowed by 
gauge invariance for the flattest direction @|. This gives 
$i ~ 10 16 GeV for M the reduced Planck mass, which 
would be too small for the GW computed here to be 
observable by currently proposed experiments. However, 
the fact that this term is allowed by gauge invariance does 
clearly not imply that it is indeed present, e.g. it can be 
forbidden by other symmetries. Higher-order terms lead 
to larger values of $i and therefore to GW that can be 
observable. Indeed, initial VEVs can be as high as the 



Planck scale [||, e.g. in no-scale supergravity models 0]. 

To summarize, we performed the first non-linear study 
of the non-perturbative decay of supersymmetric flat di- 
rections with lattice simulations. This explosive process 
can have important consequences in cosmology, in par- 
ticular for baryogenesis and reheating after inflation. We 
showed that it generates a stochastic background of GW 
that is very sensitive to the initial VEV of flat direc- 
tions when they start to oscillate, to their soft SUSY- 
breaking mass and to the time when the universe be- 
comes radiation-dominated. The specific spectral prop- 
erties that we computed can be used to distinguish these 
GW from other possible cosmological backgrounds. Of 
particular interest is that they fall in the Hz-kHz fre- 
quency range, where the signal is not expected to be sig- 
nificantly screened by stochastic backgrounds from as- 
trophysical sources. We showed that, even in a simple 
model, these GW can be observable by upcoming experi- 
ments, including Advanced LIGO, if some flat directions 
acquire a large enough initial VEV. This typically re- 
quires the flatness of the potential to be protected up 
to high order by symmetries, which is model-dependent. 
Our results strongly motivate the study of more complex 
models, with multiple flat directions and local gauge in- 
variance [t| [l(J. In particular, vector fields are likely to 
increase the range of parameters leading to an observ- 
able signal (la |. Any possibility of direct detection has 
to be studied in depth, because it would provide crucial 
informations on both supersymmetry and the very early 
universe. 
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